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ON TRAVELING WAVES IN BEAMS 1 

By Robert W. Leonard and Bernard Btjdiansky 


SUMMARY 

The basic equations of Timoshenko for the motion of vibrating 
nonuniform beams , which allow for effects of transverse shear 
deformation and rotary inertia, are presented in several forms, 
including one in which the equations are written in the directions 
of the characteristics . The propagation of discontinuities in 
moment and shear, as governed by these equations, is discussed. 

Numerical traveling-wave solutions are obtained for some 
elementary problems of finite uniform beams for which the prop- 
agation • velocities of bending and shear discontinuities are 
taken to be equal . These solutions are compared with modal 
solutions of Timoshenko } $ equations and, in some cases, with 
exact closed solutions. 

INTRODUCTION 

The theoretical analysis of transient stresses in aircraft 
wings and fuselages subjected to impact loadings has gen- 
erally been performed by means of a mode-superposition 
method that uses the natural modes of vibration predicted 
by the elementary engineering theory of beam bending. 
(See, for example, refs. 1 to 3.) For very sharp impact 
loadings, however, this approach is known to have certain 
shortcomings. For sharp impacts of short duration, many 
modes are often required for a satisfactory degree of con- 
vergence (see, for example, ref. 4); in addition, the use of 
elementary beam theory in the calculation of the higher 
modes of vibration is inaccurate because of the neglect of, 
among other factors, the effects of transverse shear deforma- 
tion and rotary inertia which become increasingly important 
for higher and higher modes (ref. 5). 

A classically recognized alternative to the modal method 
of calculating transient stresses in elastic bodies is the travel- 
ing-wave approach, which seeks to trace directly the propa- 
gation of stresses through the body (ref. 3). Although the 
traveling-wave concept has been successfully used to treat 
such simple problems as longitudinal and torsional impact 
of rods, only recently have serious attempts been made to 
study the transient bending response of beams by this 
approach. Fltigge (ref. 6) was apparently the first to point 
out that elementary beam theory could not serve as an 
adequate basis for the traveling-wave approach since the 
elementary theory predicts that disturbances propagate with 
infini te velocity; he showed, however, that a traveling-wave 
theory could be constructed by modifying the elementary 


theory, as Timoshenko (ref. 7) did, to include first-order 
effects of transverse shear deformation and rotary inertia. 
On the basis of this more accurate theory, Fltigge found that 
discontinuities in moment and shear travel along the beam 
with finite, and generally distinct, velocities. A similar 
analysis was carried out independently by Robinson (ref. 8) 
who, exploiting the hyperbolic nature of Timoshenko's 
equations, proposed the use of approximate methods of 
solution and gave some numerical results for a particular 
example. Pfeiffer (ref. 9) also suggested the possibility of 
step-by-step solutions by the method of characteristics. In 
reference 10, Uflyand attempted an analytical solution of 
Timoshenko's equations for the case of a simply supported 
beam subjected to a sudden application of load; however, as 
was shown by Dengler and Goland (ref. 11), Uflyand's work 
is marred by the fact that he applied boundary conditions 
that are incorrect for Timoshenko's theory. The only 
known example of an exact traveling-wave solution based on 
Timoshenko's theory was presented by Dengler and Goland 
fpr the case of an infinitely long beam subjected to a con- 
centrated impulse. 

Thus, although the use of Timoshenko's theory as a basis 
for the transient-stress analysis of beams has been seriously 
considered, few problems have actually been solved. Addi- 
tional basic studies of Timoshenko's equations and their 
solution, particularly for finite-length beams, constitute 
necessary prerequisites to the development of practical 
methods of dynamic-stress analysis based on the traveling- 
wave approach. In the present report, several ■ specific 
problems of transient loading of uniform beams of finite 
length are considered and their solutions by various proce- 
dures, all based on the Timoshenko theory, are presented. 
These procedures are (a) numerical step-by-step integration — 
the “method of characteristics," (b) mode superposition, and 
(c) exact closed-form solution. The examples are all for the 
special case of equal propagation velocities of shear and 
bending disturbances; only for this case have exact solutions 
been found in closed form. For the sake of completeness, 
the presentation of these solutions is preceded by an exposi- 
tion of the basic equations of Timoshenko's theory, a deriva- 
tion of the characteristic lines and characteristic forms of 
these equations, and a discussion of their implications con- 
cerning propagation of disturbances. 


i Supersedes NAOA TN 2874, “On Traveling Waves in Beams 1 ' by Robert W. Leonard and Bernard Budlansky, 1953. 
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cross-sectional area 
Young’s modulus of elasticity 
shear modulus of elasticity 
cross-sectional moment of inertia 
length of finite beam (arbitrary length for infinite 
beam) 

internal bending moment (see fig. 1) 
dimensionless internal bending moment, ML/EI B 

-T 

TiJ 


vertical shear force on a cross section (see fig. 1) 
dimensionless vertical shear force, VL*/EI B 
propagation velocity of bending discontinuities, 


V 

propagation velocity of shear discontinuities,-^ 


EIb 

pit 


pA t 


operator used in the Laplace transformation 
intensity of distributed external loading 
dimensionless intensity of external loading, gL 3 jEI B 
cross-sectional radius of gyration 
time 

velocity of deflection, y t 
dimensionless velocity of deflection, v/cn 
coordinate along beam 
deflection (see fig. 1) 
dimensionless deflection, y/L 

denotes a discontinuity in quantity immediately 
following 


dimensionless coordinate along beam, xjL 
density of beam material 
dimensionless time, CitjL 

rotation of cross section about neutral axis (see 

1) 

velocity of rotation of cross section, \[/ t 
dimensionless velocity of rotation, coL/ci 


Subscripts: 

B contributing to resistance of beam to bending 

i contributing to inertia 

S contributing to resistance of beam to shearing 

Xjt&r indicate partial derivatives with respect to those 
quantities 



BASIC EQUATIONS 

NONUNTFORM BEAMS 

Timoshenko’s equations. — The Timoshenko theory of 
beam bending (ref. 7) constitutes a modification of elemen- 
tary beam theory that attempts to account for the effects 
of transverse shear deformation and rotary inertia; the basic 
assumption of elementary theory — that plane sections 
remain plane — is retained. The moment M, shear V, de- 
flection y, and cross-sectional rotation t p (see fig. 1) of a 
nonuniform beam subject to a dynamic lateral loading q 
are governed, according to this theory, by the following 
four simultaneous partial-differential equations: 


A£-\-EI B ip s — 0 (la) 

V-AsQty -0 = 0 (lb) 

M x V-\-pIi'Ptt = 0 (lc) 

V x —pA<y tt +q=0 (Id) 


The first two equations constitute elastic laws relating the 
deformations to the internal loading. Equation (la) ex- 
presses the same relationship between moment and cross- 
sectional rotation as that given by elementary beam theory 
Equation (lb) stipulates a linear relationship between the 
shear V and the shear angle yx—^P at the neutral axis; A a is 
the so-called “effective” shear-carrying area, different from 
the total area A { since the true shear angle actually varies 
over the cross section. Equations (lc) and (Id) prescribe 
rotational and translational equilibrium, respectively, with 
the term plpp tt representing the contribution of rotary inertia. 

The moment and shear may be eliminated from equations 
(1) to yield two simultaneous partial-differential equations 
in y and \p: 


(EI B \p^) x +A 8 G (y x —ip)— pltftt^ 0y 
[A 3 G (y*— 4d]x— pAiy tt + 2=0 ) 


( 2 ) 


This form is convenient for finding the normal modes and 
frequencies of free vibration (2=0) predicted by Timo- 
shenko’s theory and for carrying out modal analyses that 
make use of these modes. In this theory, a natural mode 
is described by a pair of functions [y(x)j ip(x)] rather than a 
single function y(x), as is the case in elementary beam theory. 

For a traveling-wave analysis, however, it is advantageous 
to return to the original system of equations (eqs. (1) ) but 
to put them in a more convenient form by differentiating 





V+V x dx 


Figure 1. — Positive distortions and positive internal forces and moments associated with a typical beam element. 
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equations (la) and (lb) with respect to time, 
tions become: 



M t = 0 


The equa- 
(3a) 


Cl ~CT y, -“ = ° (3b) 

M x -\- pli&t — V =0 ( 3 c) 

V r — pA iVt + 2 = 0 ( 3 d) 


where the new variables, linear velocity v and angular 
velocity co, have been introduced for y t and \f/ t) respectively. 

Equations (3) comprise four first-order linear partial- 
differential equations in the four variables v , co, M, and V. 
Furthermore, equations (3a) and (3c) contain derivatives 
of only M and co, whereas equations (3b) and (3d) contain 
derivatives of only v and V. These facts are exploited in 
the next section in seeking characteristic lines and char- 
acteristic forms of these equations. 

Characteristics and the characteristic form of the equa- 
tions. — In equations (3a) and (3c), the variables M and co 
are differentiated with respect to both space and time; it 
would be advantageous to replace them by equivalent equa- 
tions each involving only total derivatives (or differentials) 
in a particular direction in the space-time plane. The lines 
in the space-time plane having these particular directions — 
the characteristic lines or so-called “characteristics” — and 
the equivalent equations written in these directions are found 
as follows (ref. 12). * 

A linear combination of equations (3a) and (3c), 


) + M x +p/iCo t — V=0 (4) 


is formed, where the function p is to be determined in such a 
way that the partial derivatives in equation (4) combine to 

give total derivatives and ^ in the direction of an as yet 


unknown characteristic line [ 2 ( 0 ), £(o)]. In order that the 
terms involving derivatives of M combine in the form 




dt_dM 

do do 


the function p must satisfy the following equation: 

dt 

do_ dt 

^~dx dx 
do 

where ^ is the required slope of the characteristic line. 

Similarly, in order that the terms involving derivatives of co 
combine in the form 

dx dt_do3 
^d^'d* d* 


the function p must also satisfy the equation 
* 

Ph 


dt 

do dt 


pEI B dx dx 
do 


Since the characteristic slope must be the same in both cases,. 
p is defined by 

M ei b 

Thus, the two values of p and the corresponding characteristic 

, dt 
slopes ^ are: 

1 1 

<&_ 1 
dx Ci J 


(5> 




dt 


( 6 ) 


dx Ci J 

[ei~ 

where Then, multiplying equation (4) by dt 

and using equations (5) yields 


when 


— dM-\- pli d<o — V dt=0 

Ci 


dt=—dx 

Ci 


Similarly, using equation (6) gives 


when 


—— dM+pIt dw-V dt = 0 

Ci 


dt= dx 

Ci 


In an analogous fashion, it can be shown, from equations 
(3b) and (3d), that 


- dV-pAi (fo+(Vf-^ dt= 0 

Ca \ C2 / 

when 

dt== — dx 

C2 

and 

~ dV-pAi dv+(z -^- «) dt = 0 

when 


where C 2 




Then, the system of equations (3) has associated with it 
four real characteristic directions and is thus “totally 
hyperbolic” (ref. 12). A network of characteristics can 
readily be constructed without prior knowledge of the 
unknowns M, V , co, and v since their slopes depend only on 
the material and geometrical properties of the beam. For 
uniform beams, as well as for tapered beams having uniform 
material properties and geometrically similar cross sections, 
Ci and c* are constant, and the characteristics will therefore 
be straight; in general, however, the characteristics may be 
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Figure 2. — The characteristics of Timoshenko’s equations for a point 

in the x t t plane. 


curved. Figure 2 illustrates the four characteristics passing 
through a point P in the spacer time plane with the character- 
istics designated as 


I+: 

dd_ 1 
dx Ci 

I—: 

dt 

dx 

n+: 

dt_l 

dx Ct 

n-: 

dt_ 

dx i 


It is known that, by virtue of the totally hyperbolic character 
of the basic equations, the values of the unknowns M, V , 
a?, and v at the point P depend only on their initial values at 
£=0 between the points Xi and x$ on the beam (ref. 12). 
Furthermore, these values at P can, in turn, have influence 
only on points lying in the region above P enveloped by the 
1+ and I— characteristics through P. Thus no signal 
can proceed along the beam with a velocity greater than 
Ci (which is generally larger than c*, as illustrated in fig. 2). 

For the sake of easy reference, the four characteristic 
differential forms of the basic equations are grouped below. 

Along I+: du — V dt = 0 (7a) 

Ci 

Along I — : — dM — p/* dw-\~V dt = 0 (7b) 

Ci 

Al ong 11+ : — dV — pAi dv-\-(pA i c*o)-{‘q)dt=0 (7c) 

Along HE — : — dV-\-pAi dv-\-(pAiCiCo — q)dt= 0 (7d) 

Ci 

Related forms of these characteristic equations have been 
written by Robinson (ref. 8) and Pfeiffer (ref. 9) . 

Propagation of discontinuities. — Characteristics are lines 
across which discontinuities may exist (ref. 12) ; indeed, this 
property is often taken as the basic definition of a character- 
istic. In the present problem, discontinuities (or jumps) in 


M and co can therefore exist across the 1+ and I — char- 
acteristics, and discontinuities in V and v can exist across 
the H+ and II — characteristics. Hence, a jump in M 
or co will propagate with the velocity Ci, whereas a jump in V 
or v must proceed with the velocity c%. It should be remarked 
that such discontinuities would appear only in the limiting 
case of a beam subjected to an instantaneous loading. The 
solution of such idealized problems, which are often instruc- 
tive, requires a knowledge of the laws governing the varia- 
tions in strength of these discontinuities as they propagate 
through the beam. These laws are determined below for 
nonuniform beams for which it is assumed that the condition 
c 1 —csx does not hold over any portion of the beam; in other 
words, the I and H characteristics are distinct. The 
special case where ci=c^ is considered subsequently when 
uniform beams are discussed. 

Let a and b be two points on a I— characteristic on either 
side of a particular 1+ characteristic 



If M is discontinuous across the 1+ characteristic, then 
M a — M h retains a finite value 8M as a and b are allowed to 
approach the 1+ characteristic from either side. Con- 
sequently, from equation (7b) written along the I — 
characteristic, 

-SM—pI ( Sw=0 

Ci 


since dt approaches zero as a and b approach each other. 
Thus, everywhere along a 1+ characteristic, jumps 8M and 
5co across this characteristic are related by 



8M= Cipli 8a> 

(8) 

Similarly, the jumps across the other characteristics can be 
readily shown to satisfy 

Along I— : 

8M= — Cipli 5 co 

(9) 

Along II+: 

8V— — &ipAi8v 

(10) 

Along II— : 

- 8V=c^pAi8v 

(ID 


A jump in M is thus always accompanied by a definite jump 
in co; similarly, jumps in V and v are always coupled together. 

The variations in the magnitude of a discontinuity as one 
proceeds along a characteristic may be determined in the 
following manner. Equation (7a) is written for the upper 
side and the lower side of the 1+ characteristic; then, since 
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V is continuous across a 1+ characteristic, the difference of 
the two equations yields 

d(SM) -\-C\pIid{Su) = 0 

along I+. Eliminating Sw by using equation (8) gives 
d(8M)+CipIid (£r)=0 

By carrying out the indicated differentiation in the second 
term and dividing by 2 8M, the following result is obtained: 

d(8M) 1 djciplj ) 

8l \1 2 Ci pli 

Solution of this equation gives 

as the relationship between the magnitude of the jumps in 
M at two points 1 and 2 on the 1+ characteristic. 

It can be shown that the identical relationship holds be- 
tween jumps in M at two points on a I— characteristic. 
Similarly, it can be found that, on 11+ and II— character- 
istics, 

(I3) 

for any two points 1 and 2 on a given characteristic. The 
corresponding variations of the jumps $a> and 8v are, of 
course, readily determined from equations (8) to (11). 


UNIFORM BEAMS 

Nondimensional form of the equations. — The examples to 

be presented in this report are all concerned with beams 

having uniform cross-sectional size and shape and uniform 

material properties throughout their length. For such 

beams, it is convenient to express Timoshenko's equations in 

nondimensional form. Thus, equations (3) may be written 

• a t t? ML VL S - coL -j - v 
m terms of V-nT"* «= — > and as 

iLIb -CjIb Ci C? 

5 { +M r = 0 (14a) 

(14b) 


M(+Zi T -V= 0 (14c) 

Vt-R - F t + 2=0 (14d) 

C2 


where r=^, and The quan- 

tity L refers to the beam length for all beams except 
those of infinite length, in which case any convenient arbi- 
trary length may be chosen for L. 

The characteristics of Timoshenko's equations for a uni- 
form beam are defined in the £,t plane by the families of 
straight lines 


I+: 


dr 

d£ 


1 


I—: 

^=-1 

di 

(15b) 

n+: 

dr Ci 

d£ c. 

(15c) 

n-: 

<3jS 

II 

(15d) 

The nondimensionalized characteristic forms 
equations become 

of the basic 

Along I+: 

dM-\-dto — V dr= 0 

(16a) 

Along I— : 

dM — cK+F dr= 0 

(16b) 

Along 11+ : 

dV-R <fip+(.R5+| 2 ) dr= 0 

(16c) 

Along II — : 

JV+R dF+(tf5-| 2 ) dr =0 • 

(16d) 


In addition to the restriction to uniform beams, for which 
Ci and Ci are constants, the examples presented herein are 
further limited to those beams for which the propagation 
velocities Ci and c* are equal. This assumption has been 
made because the simplifications that result not only permit 
the ready attainment of numerical solutions but also, in par- 
ticular cases, permit the attainment of exact closed solutions 
for comparison. Since, for this very special case, the char- 
acteristics II coincide with the characteristics I, equations 
(16) may now be written, fof #=0, as 


Along I+: 


Along I— : 

where A=-^ 
2r< 


dM-\-dd> — V dr=0 
dV—4\ 2 dv +4A 2 w 0 


dM-dL+Vdr= 0 
dV +4\ 2 <fr +4A 2 a> dr= 0, 


(17a) 

(17b) 


Propagation of discontinuities when ci=C 2 . — Equations (8) 
to (13), which describe the behavior of discontinuities in a 
nonuniform beam, also describe, as a special case, the be- 
havior of discontinuities in a uniform beam for which 
They show that such discontinuities propagate with constant 
magnitude. 

However, when the beam has properties such that Ci=c*, 
these equations are no longer valid. Equations which are 
valid in this case may be derived in precisely the same way 
by using equations (17) instead of equations (7). Discon- 
tinuities in such beams can be shown to be related by 


Along I+: 


(8M=8co 
UF=— 4 \ 2 $F 



(18a) 


(16a) 


Along I— : 


(18b) 
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as they propagate, and they can be shown to vary in mag- 
nitude according to the equations 

Along 1+ and I—: < / (19) 

( 5 ^+ 2x2 ^=°) 

Equations (19) may be solved simultaneously to obtain 
( 8 00 = A COS (Xr — JB ) ) 

Along 1+ and I— :< _ > (20) 

(bV=—2\A sin (Xr — B) ) 

where A and B are arbitrary constants which must be eval- 
uated by using known values of 5co and dV at some point on 
the characteristic. The variations in bM and bv can then 
be readily found by using equations (18). 

Thus, for the case Ci=C 2 , discontinuities in a uniform beam 
do not propagate unchanged but vary in magnitude sinu- 
soidally as they progress through the beam. 

LIMITATIONS OF THE THEORY 


practical structures; the last procedure, although exact, 
would rarely be useful in practice and is introduced herein 
chiefly as a check on the accuracy of the first Wo. 

The detailed descriptions and applications of the throe 
methods are contained in appendixes A, B, and 0. In 
brief, the numerical procedure exploits a grid of character- 
istic lines as shown in figure 3 (a). For the case ci—c? that 
is under consideration, this grid consists of Wo families of 
lines in the £,r plane with slopes 1 and —1. A recurrence 
formula^is developed in appendix A that gives the values of 
w and V at station 1 of a typical interior mesh (see fig. 3 (b)) 
in terms of the values of a> and V at stations 2, 3, and 4. 
Repeated application of this formula, together with the 
use of special formulas for the half-meshes at either end of 
the beam _and the knowledge of the magnitudes of jumps 
in co and V across characteristics where they occur, permits 
co and V to be calculatedjhroughout the £,r plane. Subse- 
quent determination of M and v is achieved by means of 
simple addition formulas utilizing these calculated values 
of co and V. 

Although the solutions derived in appendix B have actually 
been obtained by Laplace transform techniques, they have 


It may be well to insert a word of caution about the 
applicability of Timoshenko's theory. The investigations of 
Prescott (ref. 13) and Cooper (ref. 14) have shown that, 
when the response of a beam includes components whose 
wave length is small compared to the depth of the beam, 
the assumption of Timoshenko's theoiy that plane sections 
remain plane after 'bending is, as might be expected, too 
restrictive. Since applied disturbances which could give 
rise to discontinuities would obviously excite even the short- 
est wave length in the spectrum of the response, the results 
obtained by application of Timoshenko's theory to such 
hypothetical problems cannot, in themselves, have practical 
significance. However, carrying out solutions involving dis- 
continuities is a useful means of testing methods of solution 
of the Timoshenko equations with a view to applying these 
methods to problems in which discontinuities do not exist. 
Furthermore, the admittedly inaccurate response to an in- 
finitely abrupt disturbance may be used to obtain the correct 
response to disturbances of a more practical nature through 
the application of Duhamel's superposition integral. 

SPECIFIC EXAMPLES — FINITE UNIFORM BEAMS WITH ci=C2 

Three specific examples are considered; they are: a canti- 
lever beam given a step velocity at the root, a simply sup- 
ported beam subjected to a step moment at one end, and 
a simply supported beam subjected to a ramp-platform 
moment at one end. 

METHODS OF SOLUTION 

In the examples, the results of calculation by the following 
three methods are compared: (a) numerical step-by-step 
integration along the characteristics, (b) normal-mode super- 
position, and (c) exact closed-form solution. The first Wo 
procedures are approximate in character, but they could 
conceivably be generalized sufficiently to be applied to 




(a) Space-time plane., (b) Typical interior mesh. 

Figure 3. — Characteristic grid for application of numerical procedure 
to uniform beams with disturbance applied at £=»0. cj=C|. 
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been termed “modal solutions” because they are exactly 
those which would result from an application of the usual 
mode-superposition process. The exact closed solutions in 
appendix C have also been obtained through the use of 
Laplace transforms. 

CANTILEVER BEAM GIVEN A STEP VELOCITY AT THE ROOT 

The first example to be considered, the response of a 
uniform cantilever beam given a step velocity v=l at 
the root, is the .equivalent of the most severe idealized 
landing problem, the instantaneous arrest of the root of a 
moving cantilever beam. Computed results for the shear 
and moment at the root of such a beam having a ratio of 
length to radius of gyration of 10 (\=5) and properties 
such that Ci=Ca are presented in figures 4 (a) and 4 (b), 
respectively. Two separate curves obtained by the numeric 
cal procedure are shown — one from a grid that divides the 
beam into 10 segments and the other from a 20-segment 
solution. The modal solution includes the contributions of 
the first eight modes. Results given by an exact closed 
solution are shown for both the shear and moment at the 
root up to the time t— 2. These exact results are actually 
those for an infinitely long beam, since the influence of the 
free end is not felt at the root until r=2. After r=2, the 
influence of the free end is felt and, in this case, an exact 
solution is not feasible. To illustrate the time range covered 
in the plots, the point corresponding to half the period of the 
first mode of vibration of the beam is indicated on the time 
scale of each plot. 

In figure 4 (a), the shear discontinuities evident in the 
numerical solutions occur each time the discontinuous wave 
front returns to the root after being reflected at the free end. 
The beam is seen to react violently to each of these boosts by 
the wave front, with more and more oscillations occurring 
after each succeeding jump. The frequency of these oscilla- 
tions tends to increase with each succeeding boost until 
limited by the finite time interval. In these regions of violent 
oscillation the accuracy of the numerical solutions is ob- 
viously questionable; indeed the question arises as to whether 
these oscillations are really predicted by the theory or are 
merely the result of some instability in the numerical process. 
This question is resolved in the next section in which the 
simply supported beam is considered. At any rate, away 
from the regions of violent oscillations, the comparisons with 
the modal solution are favorable, and, for r<2, the fine-grid 
numerical solution almost coincides with the exact closed 
solution valid in this region. 

The comparisons between the numerical and modal solu- 
tions are very good in figure 4 (b) where the time history of 
the moment at the root is plotted. Again, the fine-grid nu- 
merical solution nearly coincides with the exact closed solu- 
tion for r<2. 

SIMPLY SUPPORTED BEAM WITH AN APPLIED END MOMENT 

Step moment. — A simply supported unif orm beam with a 
ratio of length to radius of gyration of 10 (X=5) and prop- 
erties such that Ci=C 3 is subjected to a step moment M= 1 


at the end £=0. Computed results for the shear at $=0 are 
presented in figure 5 (a) and the time history of moment at 
the center of the beam is presented in figure 5 (b). The 
point corresponding to the full period of the first mode is 
marked on the time scale of each plot. The numerical curves 
for both shear and moment were obtained by a 20-segment 
solution. The modal curves were obtained by adding dy- 
namic corrections to the static solutions, the dynamic correc- 
tions being calculated with six modes for both the shear and 
the moment. 

This example affords an answer to the question raised in 
the preceding section with regard to the stability of the nu- 
merical procedure following the passage of a discontinuous 
wave front. The fact is that the violent oscillations that 
occur after the discontinuity actually appear in the exact 
solution (fig. 5 (a)) and are hence inherently present in the 
theory. 

In figure 5 (a) the scale of dimensionless vertical shear 
happens to be precisely the dynamic overshoot factor; it is 
of interest to note that values at least 15 times the static 
shear are predicted when the moment is applied suddenly. 

For shear (fig. 5 (a)), the inaccuracies in the numerical 
solution just after the discontinuities are evident; however, 
the numerical results approximate the exact solution very 
well elsewhere. A similar observation may be made for the 
modal solution, which, as would be expected, ignores the dis- 
continuities and violent oscillations caused by the wave front. 
The numerical and modal curves for the moment in figure 
5 (b) follow the same pattern, agreeing well everywhere ex- 
cept in the regions immediately following the discontinuities. 

Ramp -platform moment. — Computed results for the sh ear 
at the end £=0 and the moment at- the center of the same 
simply supported uniform beam subjected to the applied 
ramp-platform moment 

A?(0,t) = t (O^T^l) 

•2Z(o, t) — i (r>i) 

are presented in figures 6 (a) and 6 (b). Again, the period 
of the first mode of vibration is marked on the time scales. 

The numerical curves were recomputed with a 20-segment 
grid for the new forcing function; however, the results for the 
modal and exact solutions were obtained by means of a 
superposition of the preceding results. This superposition 
was carried out analytically for the complete modal solution 
and for the exact solution in the range t< 2. In the range 
t> 2, it was necessary to carry out the superposition for the 
exact solution numerically. 

In figure 6 the time to peak value of the applied moment 
is seen to be approximately one-seventh the period of the 
first mode; predicted shear values approximately three and 
one-half times the static response occur. 

With the removal of the discontinuity, there remain no 
high-frequency oscillations which the numerical solution 
might be unable to represent. In fact, all three solutions 
for the shear and both moment solutions are seen to be in 
excellent agreement. 
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(a) Time history of the non dim ensional vertical shear at the root. 

(b) Time his tory of the nondimensional moment at the root. 

Figurb 4. — Response of a uniform cantilever beam subjected to a step velocity at the root. 
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Figure 6. — Response of a uniform simply supported beam subjected to a ramp-platform moment at £=>0. 


CONCLUDING REMARKS 

Timoshenko’s equations for the motion of vibrating non- 
uniform beams may be written in a characteristic form which 
appears to be well-suited to solution by numerical methods. 
In the examples presented in this report, all of which are for 
uniform beams with equal propagation velocities of bending 
and shear discontinuities (^1=02), the solutions by the numer- 
ical and modal methods generally agree well with each other 
as well as with exact closed solutions where these have been 
obtained. However, the modal method, of course, fails 
entirely to reproduce discontinuities in shear or moment, 
and the numerical procedure, although it yields these 
discontinuities, represents the initial oscillations which 
follow with decreasing accuracy as more and more reflections 
of the wave front occur. In the more practical situations 


where discontinuities do not exist, these difficulties will, of 
course, not arise. 

The results of the examples carried out by the numerical 
traveling-wave procedure encourage the viewpoint that 
such traveling-wave analyses may eventually be of practical 
usefulness. This kind of procedure is inherently sinlple and 
straightforward and has the advantage that the bulk of the 
labor involved is routine computation rather than mathe- 
matical analysis. It should be emphasized, however, that 
numerical solutions of Timoshenko’s equations have been 
demonstrated only for uniform beams in which the propaga- 
tion velocities c x and <h are equal; numerical procedures for 
the general case where they are unequal remain to be devel- 
oped and tested. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 28 , 1952 . 



APPENDIX A 

NUMERICAL SOLUTIONS FOR UNIFORM BEAMS WITH ci=C2 


MATRIX FORMULATION 

Let the differential equations (17) be replaced by the 
finite-difference equations 


Along I+: 


Along I— : 


AM+Aco — FAt— 0 

(Ala) 

AF- 4X*A© +4X*5At= 0 

(Alb) 

AM— Aw-(- FAt= 0 

(Ale) 

AF+4X*A© +4X 2 «Ar= 0 

(Aid) 


which may, in turn, be subtracted to obtain the single 
equation 


b>i— co*— CO3+CO4 — ~{Vi— F<) = 0 (A10) 

Similarly, equations (A8) and (A9) may be subtracted from 
equations (A4) and (A5), respectively, to obtain 

Vi~ 2 P 2 + 4\ 2 (th— V 4 ) +2X a Ar(w x — co*) = 0 

and 


and consider a closely spaced network of 1+ and I— char- 
acteristics in the space-time plane as shown in figure 3 (a). 
Let the corners of a typical interior mesh of this grid be 
designated as shown in figure 3 (b). A step-by-step integra- 
tion formula for co and T r may now be derived in the following 
manner. 

Equations (Al) may be written along the upper sides of 
the typical mesh to obtain 

M-M 3 +5i-5 s -y(F 1 +F 2 ) = 0 (A2) 

Ht-Mz-Zi i+S3+y(Fi+F 3 ) = 0 (A3) 

Vi— Vi— 4X s (Fi— p s )+2X 3 Ar(Si+w 2 ) = 0 (A4) 

Vi~ F 3 +4X*(©i- v 3 ) +2\ 1 At(Z 1 +Z 3 ) = 0 (A5) 


Vi — 2 F 3 + F 4 -(- 4 X*(©i — ©4) + 2 X s Ar(«i — W4) — 0 
and these may be added to derive 


V 1 - Vs- V 3 + F 4 +2X i Ar(^ l - «4) = 0 (Al 1 ) 


Equations (A10) and (All) may now be solved simultane- 
ously to obtain the step-by-step integration formula 


«1 


1 1 

(X Ar)*+1 (j/j Vs+V 3 



where 

1 — " 

[A»]= 2 

2 X*Ar 1 _ 

and 


(A12) 


where « and V have been assumed to vary linearly in the 
small intervals between the comers. Obviously equations 
(A2) to (A5) may be solved simultaneously to obtain the 
four quantities V, u, v, and M at point 1 in terms of their 
values at points 2 and 3. However, it is noted that, if F, 
w, r, and A? at points 2 and 3 were determined by a similar 
process for the preceding meshes, they already satisfy the 
equations 


(A6) 

(A7) 

(AS) 
(A9) 

Equations (A6) and (A7) may be added to equations (A2) 
and (A3), respectively, to obtain 

'M\ — 'Mi-\-<d 1 — 2co 3 -|-^4 — 2*(F[ F4) = 0 

and 

Mi— Mi— aJi+2w 3 — £ 4 -|-y (Fi- Fi) = 0 


M- Mi- Zs+ Zi+^-{Vs+ F4) = 0 

M 3 -Mi+Zs-ui~(Vs+Vi)=0 

Vs- Vi+ 4 \ 3 (v- 2 - ©4) +2X s Ar ( 5 a + W4) = 0 
F 3 — Vi— 4 X*(© 3 — ©4) +2X i Ar (aj 3 + W4) = 0 



(X At)* — 1 
4X s Ar 


— Ar 
(X At)* 


J 


Thus co and V may be determined at every interior mesh 
point in the space-time plane by the repeated application 
of formula (A12) to each mesh as it is encountered. The 
half-meshes which occur at the vertical left and right bound- 
aries of the plane (fig. 3 (a)) require special formulas incor- 
porating the known boundary data. These formulas may 
be derived from equations (Al) by a procedure similar to that 
used in obtaining equation (A12). Boundary formulas for 
some specific examples are presented in subsequent sections 
of this appendix. 

Besides boundary data, initial data must be provided in 
order that a step-by-step solution may be begun. In all the 
examples considered, the beam is initially at rest and has a 
disturbance applied at the point £=0 beginning at time 
r=0. The region r<£ (without grid lines in fig. 3 (a)) 
therefore is one of zero stress and motion and co and *V are 
given along the line r=£ by the known conditions at the 
wave front. It is with these values that each numerical 
solution is begun. 

Discontinuities in co and V offer no special difficulties since 
they propagate so as to be always located on characteristics 
which, of course, can be made part of the basic grid. Thus 
they are simply added as they are encountered. 
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Once 3 and V have been determined at the grid intersec- 
% tions by the step-by-step process, A? and v may be found by 
repeated application of equation (A2) or (A3) and equation 
(A4) or (A5) along the proper grid lines from boundaries 
where TI and v are known. For example, if A? is known 
along the left boundary, it may be found at a point g by 



Typical left-boundary and right-boundary half-meshes are 


^(0 , r )= 1 1 
*(0, r) = OJ 


Apphcation of equations* (Alb) and (Aid) to the sides of 
the left-boundary half-mesh and proper combination of the 
resulting equations to eliminate v produces 

V=2F 3 -y, (A16) 

Similarly, apphcation of equations (Ala) and (Ale) to a 

right-boundary half-mesh produces 

CO! =2w 2 — 3* (A 16) 

so that a complete set of integration formulas for determin- 

ing co and V for this problem is now available. In addition, 
if the unused equations (eqs. (Ala) and (Ale) for the left- 
boundary half-mesh and eqs. (Alb) and (Aid) for the right) 
are applied and combined to eliminate A? and v at points 3 
and 2, the recurrence relations 

(Fi-V) (A17) 



At 

1 



l7(l,r)=0 


applying equation (A2) successively to the intervals ab, be, 
cd, . . . fg. The resulting expression for M t becomes 

M t =M a +5 a -s I +Ar(Iy a +y 6 +. . .+yrt-|y,) (Ai3) 

This procedure is seen to correspond to integration of 
the first of differential equations (17a) by means of the 
trapezoidal rule. 

SPECIFIC PROBLEMS 

Cantilever beam given a step velocity at the root. — If the 
root {=0 of a uniform cantilever beam is given a step 
velocity t> = l at time r=0, the boundary conditions may be 
written 

v(0,r) = l\ 

3 (0,r) — 0 ( 

_ } (Al4) 

Af(l,r) — 0 i 

F(1,t) = o) 


. for left-boundary half-meshes and 

2\a 1^2 d rjr Gh—S*) (Al8) 

for right-boundary half-meshes are obtained. These equa- 
tions may be used to compute M at the root and v at the 
tip after 3 and V have been determined everywhere. 

Since the applied disturbance is initially discontinuous, 3 
and V will be discontinuous along the line r={ (see fig. 
3 (a)). In order to begin the step-by-step solution for 3 and 
"V, these discontinuities must be determined in advance. 
Furthermore, since the wave front is reflected back into the 
beam at either end, 3 and V* will also be discontinuous along 
the lines r=2 — £, r=2+£, . . and 63 and 6V must bo 
determined along each of these lines before the step-by-step 
solution can be extended beyond it. The discontinuities at 
the wave front are determined as follows. 

From the boundary conditions (A14), it is seen that 
63(0, 0)=0 and 5t?(0,0) = l. Thus, from equation (18a), 
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aV(0,0) = — 4X a , and equations (20) become, for the line r=f, 


««(£,£)=— 2Xsin xn 
5V(M) = -4 :\ i cos 


(A19) 


At£=l, the discontinuities across r=£ are 53 (1,1) —sM(l, 1) = 
— 2X sin X and 5V'(1,1) = — 4X 2 cos X, so that, if the boundary 
conditions M( 1 ,r) = V(l,r) =0 are to be satisfied, jumps must 
occur across r= 2— £ at the point (1,1) with the magnitudes 
53(1,1) = — 6M(l,l) =— 2X sin X and 5'L(l,l)=4X 2 cos X. In 
view of these initial conditions at r=l, equations (20) 
become, for the line r=2— f, 


beam into 10 and 20 segments, respectively. (The compu- 
tations for the 20-segment solution were made on the Bell 
Telephone Laboratories X-66744 relay computer at the 
Langley Laboratory.) The computed time histories of shear 
and moment' at the root have been plotted in figures 4 (a) 
and 4 (b) , respectively, up to r= 10. In order toobtain these 
time histories, the computations for 3 and V had to be 
carried through the first nine triangular regions. 

Simply supported beam with an applied end moment. — If 
the end £=0 of a uniform simply supported beam for which 
Ci=C 2 is subjected to an applied bending moment M 0 (t) 
beginning at time r=0, the boundary conditions are 


53(£,2— £) = — 2X sin X£ 
S7(Z, 2— £)=4X 2 cos X£ . 


(A20) 


Initial jump values 53 (0,2) =0 and 5V’(0,2)=4X 2 for the line 
r=2-f£ may be found from equation (A20) by satisfying the 
boundary conditions 3(0,r)=0 and »(0,r) = l. From these 
initial values, the discontinuities across r=2+£ are found to 
be the negative of those across r=£, or 


53(£,2+£)=2X sin X£1 
s7 (£,2 + £) = 4X 2 cos X|) 


(A21) 


Then, it must be true that 

^(£,4— £)=2Xsin X£ | 
bV (£,4— £) = — 4X 1 cos X£) 


(A22) 


and so forth, with the values on each succeeding line t=ti±£ 
repeating the values on the line'r=(7i— 4)±£. The 'varia- 
tions in the magnitudes of bu> and bV as the wave front 
propagates back and forth through the beam are thus clearly 
defined, and, since r<£ is a region of zero stress and motion, 
equations (A19) define the values of c3 and Y along the line 

With S and V known along r=£, the solution for co. and V 
may be begun by applying formula (A15) to the half-mesh 
in the lower comer of triangular region (7) (fig. 3 (a)). It is 
continued step by step throughout the triangle, with formula 
(A 12) being used for all interior meshes. When the first 
triangle is complete, the known jumps along the line r=2— £ 
may be added to the values computed for the under side of 
this line, and the solution may then be carried out in tri- 
angular region ® beginning again at the lower comer this 
time with formula (A16). In this way the solution may be 
carried through as many triangular regions as desired. 

Two sets of computations, with time intervals Ar=0.1 and 
At= 0.05, have been made for a cantilever beam for which 
X=5. These time intervals correspond to grids dividing the 


M (0, r ) = Mo (r) 
v (0,r) = 0 
M(l,r) = 0 
5(l,r) = 0 


(A23) 


Typical left-boundary and right-boundary half-meshes for 
this beam are 



Equations (Ale), (Ala), (Aid), and (Alb) may be written 
along the sides of the left-boundary half-mesh to give, 
respectively, 

(M 0 )i-M 3 - 3i + W 3+Y (Vi+V 3 ) = 0 

57*— C®5)4+«3— «4- ; y (F 3 +F4) = 0 

V \ — 4\ a t?3-)-2X a Ar(co 1 -l- 6J 3) = 0 

V 3 — F 4 — 4X^3 H-2X 3 Ar(5 3 +5^ = 0 

Adding the first two equations and subtracting the last two 
equations produces 

(Mo) I — (Mo) 4 CO!+ 2 ciJ 3 — (T 7 ! — VJ = 0 


and 

Y 1 - 2 Vz+ T ? 4 +2X 2 At(g 3 1 — co^) = 0 
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which may he solved simultaneously for gjj and Vi to obtain 


(X At) 2 +1 


1H S + [>]fe 


i 

— 2X 2 At|. 
(A24) 


where [A a ] and [Aj have been defined in conjunction with 
equation (A12). 

A similar process produces, for the right-houndary half- 
mesh, the formula 

Sl+Hll) <**> 


tOi 

Vi 


(XAr) s +l 


IH 


Now consider the case in which the applied moment is a 
step moment; that is, M 0 (t) = 1. Equation (A24) reduces to 


Wi 

Vi 


'(XAt) 2 +1 


![ A -lrj+Hlrj} <"> 


for left-boundary half-meshes. A similarity between the 
boundary formulas (A25) and (A26) and the interior formula 
(A12) is apparent; indeed, formula (A12) may be used 
everywhere if the sums, of the values at points 2 and 3 are 
simply replaced by twice the values at the inboard point of 
each boundary mesh. 

For this problem, as in the preceding one, the magnitudes 
of the discontinuities $&> and 8V at the wave front must>be 
determined in advance for every point on the grid lines 
defining its position. From the boundary conditions and 
equations (18a) it is apparent that So (0,0)= $^(0,0) =1 and 
5T / (0,0)=0. Thus, for the line r=f, equations (20) become 


So >= cos Xr 
SV=— 2X sin Xr 


(A27) 


Furthermore, it is found that, when the wave front is re- 
flected at either of the simply supported ends, the signs as 
well as the magnitudes of the discontinuities remain un- 
changed, in contrast to the behavior of a discontinuous wave 
front reflected at a free or fixed end as in the preceding case. 


Thus equation (A27) must be valid for the entire zigzag path 
defining the position of the wave front. 

With the discontinuities known, the step-by-step solution 
can be begun as before. Again, it is convenient to complete 
each traingular region before proceeding to the next. 

Computations for this problem were made on the Boll 
computer, with a 20-segment grid, for a beam with X=5. 
The quantities a> and V were found in the first eight trian- 
gular regions; in addition, the moment at the center was 
computed. The time histories of shear at the end £=0 and 


of moment at the center f=— have been plotted in figures 


5(a) and 5(b), respectively. 

Now consider the case in which a ramp-platform moment 
is applied to the beam; that is, M 0 (r) is defined by 


Mo (r) = T (O^r^l) | 

= 1 (r> 1) ) 


(A28) 


Equation (A24) reduces to 


«i. 

Vi 


(XAr)*+l 


fHKI+H 


o> 4 

^4 


+2Ar! 


1 

— 2X a Ar 


(A29) 


for the left-boundary half -meshes for which ngl and reduces 
to equation (A26) for all the rest. Equations (A12) and 
(A25) are, of course, still valid for interior meshes and right- 
boundary half-meshes, respectively. 

Since 5S(0,0)=$y(0,0)=0 for this case, the solution 
begins with initial values Zi=V=0 along the line r=£ and 
there are no discontinuities to be added. Thus the solution 
may progress in any convenient manner without regard for 
the position of the wave front. 

The same quantities were computed for this problem that 
were obtained in the preceding case. As before, it was 
assumed that X=5, and the calculations were performed on 
the Bell computer with a 20-segment grid. The resulting 

time histories of shear at {=0 and moment at are 
presented in figures 6(a) and 6(b), respectively. 



APPENDIX B 

MODAL SOLUTIONS FOE UNIFORM BEAMS 

The solutions carried out in this and the following appendix make considerable use of Laplace transform techniques; all 
the Laplace transforms used, most of which were taken directly from references 15 and 16, are given in table I for the 
sake of easy reference. 


TABLE I.— LAPLACE TRANSFORMS 


Number 

F(r) 

(a) 

/0> ) = J^“e-^F( r ) dr 

Reference 

1 

Fi (i 6 ) F 2 (0) dd 

(r>0) 

/ i(p)/i(p) 

15 

2 

F( T-& 

(r>{) 

e^/Cp) 

15 

3 

eP"F(r) 

(r>0) 

f(p— i\) 

15 

4 

i' p 

(r>0) 


15 

5 

tF(t) 

(t>0) 

-/'(p) 

15 

6 

„ , % d* 




F'(r) =1 -JF(r)] 

(t>0) 

P*/(p)— P^Ff+O)— f*-*F'(+ 0)— . . . — Fa-H+O) 

15 

7 

1 

( T >0) 

1^ 

15 




V 


8 

r 

(r>0) 

J_ 

15 




p* 

9 

Jo(\r) 

(r>0) 

. 1 

15 




Vp 1 +x 1 

10 

J j (Xt) 

(r>0) 

Vi p 

15 


x V Vp*+x»/ 

11 

Wo(XT)-t/,(Xr)] 

(r>0) 

KV E T 2i -) 

— 

12 

e- ttT t/o(Xr)+f/i(Xr)] 

(r>0) 

K’-V 3 ?) 

— 

13 

cos Xr «/o(Xt) “ bsin Xr 

(r>0) 

X(V^-V®) 

— 

14 

— y,( 2 Xt) 

t 

(r>0) 

Vp*+4X*-p 

15 






15 

(r>f) 

f 1 

Vp*+x 3 

15 

16 

Vt 1 -? 

(r>{) 

e -e._ e -£V^r* 

15 

17 

(r>Q 

Vr- r 


16 


a F(r)=0 for all values of r not specified. 
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EQUATIONS 


The equations of motion (2) may be written for a uniform 
beam vibrating in the absence of external distributed loading 
in nondimensionaJ form as * fc 


(?f Vrr =0 

tt —0 


(Bl) 


where y = jr The dimensionless bending moment "M and 
vertical shear V are given by (see eqs. (1)) 


Ai= — \[/( 
V=R(y ( —\f/) 


(B2) 

(B31 


If the beam is initially at rest (y(£,0)=y T (£,0)=^(£,0) = 
i/'rC^O) = 0), equations (Bl) may be transformed to (where 
transform 6 in table I has been used) 




(B4) 


¥x+R(Y t—ty— p^r=0 y 

r* cd 

in terms of the Laplace transforms F(£,p)= I e^ VT y(^r)dr 

r * cd ^ 

and SE'fop) =J 6~ PT 4 / i%, r)dr. The solutions of equations (B4) 

have the form 

Y{£,p)=A(p)e*) 


^(^,p)=B{p)^j 


(B5) 


Substituting these relations into equations (B4) leads to a 
biquadratic equation in m 

m^- 3 , s [l+(|) S ]m*+(|)*P> s +S)=0 (B6; 

and gives the following relationship between B and A: 


B=- 


p 2 


m 


(B7) 


Let the solutions of equation (B6) be written m=±imi 
and ±m 2 . Then 


m x - 




-[l+l 

©: 

l + V 

[H 

'cA*T 4R / Cl V 
\h) J V 2 Vfe/ 

2 

ks 

i+V 

1— 

©: 

4B/cA 2 
v 2 VSa/ 


and the general solution of equations (B4) is 
Y (£ ,p) = Ci cos m£+ C 2 sin mif + C& cosh + C A sinh (B8) 

^ )= i[ mi2+ (S co ^ m ^~ Gl sin m i£)+ 

p 2 J (C A cosh m&+C z sinh m 2 £) (B9) 

where Ci, 0%, Ck, and C A are functions of p which must bo 
chosen so that the boundary conditions are satisfied. 

SPECIFIC PROBLEMS 

Cantilever beam given a step velocity at the root. — If the 
root £=0 of a cantilever beam is given a step velocity 

t?=— p T = 1 at time r— 0, the boundary conditions may 
C2 

be written 


y(0,r)=^r 

Ci 

M t )=0 
&( 1 , t )=0 

^(l,-r)=i/«(l ,r) J 


(BIO) 


These boundary conditions transform to (transform 8, 
table I) 


*(0,2>)=0 
¥ { (1,P)=0 

*(l,2>)=F { (l,p)J 


(Bll) 


The constants Ci, 0%, C s , and Ci are determined by substi- 
tuting equations (B8) and (B9) into equations (Bll); the 
resulting expressions for Yfcp) and < &(£,p) may bo written 
in the forms 


and 


where 


Y(r Uy&P) 

^(t v \=Jl 
* ,P) Cip'Dfr) 


(B12) 

(B13) 


U T (Z,p)=(l — sin sinh nij-f-I cos mi cosh m^ cos 

^1+^ sin mi sinh m*-)- Q cos mi cosh m 2 ^ cosh m 2 £-|- 

( 4 sin mi cosh m* | - — 1 cos m t sinh m 2 ) sin m ^ — 

\Q . m? / 


( 


^ sin m x cosh cos mi sinh m?') sinh 77i 2 £ 


/ 
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U<t(£,p)=- 




( — sin TOi cosh. m 2 + 

m 2 L\ m i 


Q cos mi sinh m 2 ^(cos mi {—cosh {)— 

<2— ( 1— — 1 sinm,sinhm 2 +-icosmiC08hm 2 )sinm 1 {+ 
mA mi • Q, J 

sin mi sinh m a +Q cos mi cosh sinh m^J 

Z>(p)=2+^Q+^ cos mi cosh ms+ 
sin mi sinh m* 

\77ll Tflr 2 / 


Q- 






The inverse transforms j/({,r) and ^(.|,r) are determined 
by substitution of F({,p) and ^({,p) into the complex inver- 
sion integral (ref. 15). The singularities of the resulting 
integrands e Pr Y(£,p) and e Pr ^ (£,p) must therefore be exam- 
ined. Although the functions m t and are, in themselves, 
multiple-valued functions of p with branch points at p= 0, 

2 nfR — 

— * /»_ 

±iJR, and ± t— r-s, it follows from a consideration of 

the fundamental theorem of the uniqueness of solutions of 
ordinary linear differential equations that F(£,p) and it'd,?) 
must nevertheless be single valued. The integrands e pr Y (|,p) 
and £”*¥(£, p) thus have no singularities other than poles. 
The inverse transforms y({,r) and ^({,r) will therefore be 
taken as the sum of the residues at the poles of e PT Y d,p) 
and « Pr ® r d,p), respectively (see ref. 15). 

It can be shown that no singularities occur in the numera- 
tor functions Z7yd,p) and D*d,p); therefore, all the poles 
of F({,p) and it'd,?) must be introduced by the zeros of the 
denominator p J P(p). Consider first the equation 

D(p)=2 J r(^Q+-^j cos m t cosh rrh +(““^) sin m i sinh m. 2 =0 

(B14) 

This equation has an infinite number of mots p= ±p„ 
n— 1, 2, . . . oo, all of which are, in general, simple poles 
of both F({,p) and it'd,?). In addition, p= 0 can be shown 
to be a double pole of F(£,p) but not a pole of •® r (£,?)• 

The sum of the residues of c w F({,p) at the poles p=0 and 
p=±p n , n= 1, 2, 3, . . . oo, is 

y(t T )=—$ r+y^ — PZg^l’Ps) e yj JZ Y %YZlA e~ p ^~\\ 

y( Z’ T) Ci r + ^p.*L D'tol e + D'(-pl e JJ 


and, similarly, the sum of the residues of e pr ^(£,p) at p— ± p n > 
n~l i 2, 3, . . . oo f may be written 


yb(t r w£? 2A e ~ p n T 

mr) c.hip'LD'^ + TD'{-p^ e J 


where 


D'(p)=g[D(p)\ 


J<L (±-± )si 

m\m a \m 2 a 


— a — ~2 ) sin mi sinh m s~~ 


Bp'(t 



fcA 

KcJ 

2 ~ 

v \^ V/L m 2 2 (mi 4 +77ij z )J 




m 3 2 — m t 2 


771 j* 




L i +i^+ssJr“' sMi,n -- 


?iK)L 


ms 


* -i 


mf—m i s 


Mr 




J m x 2 L 


fsinmj cosh 77^2 


1 + 


However, it is seen that ^r(£,p)=Z7r(£,“-p), £/*({;, p) 
p)> and Z?'(p)=— D'(— p) so that these equations 

reduce to 

y(M-|[r+2g^»mhp„] CB15) 


and 


® I6> 


If, in equation (B14), p is replaced by ik, the result may be 
written 


2+| cosh a cos ^ sinh a sin /3=0 


L*Kt)‘ . 


(Bin 


where 


P=k 


-KD1 


■+f(D ! 

2 

i+r^v+j 

\C2/ V 


w_ 
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and fc=— ; Q being the circular frequency of vibration. 

Cl 

Equation (B17) is the frequency equation for a uniform 
cantilever beam (see ref. 5) and its solutions £=&„, n=l, 2, 
3, . . . oo, are the nondimensiona! natural vibration 
frequencies of such a beam. Equations (B15) and (B16) 
may be mitten in terms of these natural frequencies simply 
by replacing p H by ik n . The results are 



yfer)=^ T + 2 S-X’i.(0 sin &*rj 

(B18) 

and 

(B19) 

^K!,r)=2 sin *» r 

Cin=l * 

where 


*«-A[ 

“/ <v 1 > 

( 1 sinh cos /3 a cosh ct n 

A Pn 9» J 

) cosh 0,5+ 


^1— ^sin £» sinh a n +<f> n cos p n cosh cos 
(gcosh a„ sin 0,+A sinh. a* cos P^ sinh a„£+ 

(k sUl 008 |S - + *’ ^ “* “ A ) w ] 

X+ m {Q= — ^sinho^cosiSs+^.coshansin/S.jCcosjS,,?— 

cosh a n Q+<f> n ^( l+j 1 sin / 3 , sinh a„+ 

a n\ Ptt 

■7- cos 0* cosh a tt ) sinh ( 1 — — sin )S n sinh a*+ 

9* / \ “» 

4>n COS p n cosh o,^ sin /3»{J 


BJfe *( ^ 


q ~ (j ? — ~ 2 ) sin p n sinh oc n + 2 H?” - 

\j3n «./ iWg 


cos Pn cosb a„— j3. 


'(g« 2 +0 : 

Rk,* 

2 


h0 2 _ 9 l , ^ 

0 s )S, 2 (a, 






+ 


J W^Yl) 

ft.* L 1 «»W+ft, s )J) ^ " 


sin p n + 




(V W+P/) 

l_ol 

1 

mi 

bo 

d)1 

1 

to 

(1) 

• J 

a-nXan+Pi?) 


r ® i- ’© li 

L 1 ?.■(«.•+« J) ™ h .“ - 


cos A, 


The coefficients Ay w (f) and are seen to be the natural 

mode shapes associated with the frequency k n . As a check 
on the validity of the Laplace transform procedure these 
quantities may be shown to satisfy the differential equations 
(Bl) and the boundary conditions (BIO). 

Substituting y{£,r) and from equations (B18) and 

(B19) into equations (B2) and (B3) gives, for the moment 
and shear, 


m, r)«-2 ? S *5,(0 sin 

Ci n=l 


(B20) 


where 

pJ-f^Yk * 

-22,(0= — ^(l+jjjr sin sinl1 “*+ 


and 


where 


cosh a, cos cosh — ML — — sin f) n sinh a n + 

9* / \ a * 

<t>n cosh a„ cos 0*) cos &.£—(§* sinh «„ cos ft,+ 

4> n cosh a„ sin /3»)( sin sinh a„^J 


V(£>t)——2R ~ ]Lj - 2 r„(£) sin (B21) 


X8 ’® = A.[S( i+ l sin ^ 8inIi “’ + 

A cos /3, cosh a, ) sinh a*|+ ( 1— — sin p H sinh a,+ 

<p» / \ a n 

<#>« cos j3* cosh a*^ sin ft,£— (A sinh a, cos ft,+ 

cosh a, sin P^j (cos (S.5+A cosh 

The first eight terms in equations (B20) and (B21)_have 
been used to compute the quantities M(0,r) and V (0,r). 
These results have been plotted in the range O^r^lO in 
figures 4(a) and 4(b). 

Simply support ed_beam with an applied step end moment. — 
If a step moment M= 1 is applied to the end £=0 of a simply 
supported beam at time r= 0, the boundary conditions are 

y( 0,r) = 0 
y(l,r)=0 
^e(0,r) = -l 
f f (l,r)=0 

These equations transform to (transform 7, table I) 

r(o,p)=o 

F(l,p) = 0 


(B22) 


*t(0,p)=-£ 


(B23) 


*t(l,p)=0 

If conditions B23) are used to evaluate the functions Ci, 
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<7„ C 3 , and. C\ in equations (B8) and (B9), the quantities 
F(f,p) and SF(£,p) become 


7 ^ =: p (m7+m7j (cos mi * _cot 7711 8111 

cosh m^+coth m? sinh m^) (B24) 


*((,P) S 


p(mj*+m,*)L mi 


Wi 




(sin TOi£+ 


cot m i cosmiQ+ 


coth m 3 cosh m 2 £) 




m 2 


(sinh m 2 £— 


(B26) 


or they may be written in the forms 

o 328 ) 

where 

A^r(f,p)=sin mi sinh (cos cosh 7n*£)— 

cos Wi sinh sin Wif+sin m i cosh m? sinh 771 ^ 

Z?r(p) = p(m I g +^ 2 ) sin mi sinh 


and 


*(£.?)= 


_ Nf,(gp) 

D*(p) 


(B27) 


where 


Ntfa p)=mj ^m 2 2 — p*J sin Wj (cosh m 2 cosh m 2 £— 
sinh m 2 sinh m 3 f)— m 2 

(^ p 2 J sinh m 2 (cos mj cos m l |+sin mj sinmiQ 


D<t(p)=m 1 in 2 D 7 (j>) 

Here again, the functions F(?,p) and ^r(g,p) are single 
valued, and the inverse transforms y(g,r) and ^(f,r) of the 
quantities F(£,p) and '$'(?, p) are taken as the sums of the 
residues at the poles of e^Ffop) and e Pr \E r (£,p), respectively. 
Consider all the roots of the denominators Z?r(p) and Z?*(p). 
For m 1 = ±nir or nh= ±imr, n= 1, 2, 3, . . . ® (any of 
which are roots of both denominators), equation (B6) 
becomes 


P ‘+{ij + «v[i + (a)]}p- + (|)»v=o 

(n=l, 2, 3, . . . ®) (B28) 


The solutions of equation (B28) are p=±ia H and ±.ib n , n=l, 2, 3, . . . <», where 



and these points are seen to be the locations of simple poles of both F(£,p) and ^(g,p). 


The roots p— ± 


2VZ^ 


Cl 


1-(S 


of 


the equation mi 2 +77i 2 *=0 are zeros of the denominator, but these are also roots of the numerators iVr(f,p) and A 7 * (£,p) 
and are not poles of either F (£,p) ov'&{£,p). The root p=0, on the other hand, is a pole of both F (f,p) and ’J'fop) for, 
although both munerators also go to zero at this point, the denominators vanish more rapidly. The equation mi=0 has 
three roots: p= 0, which has already been discussed, and p=±i-jR. The latter points are not poles of F (£,p) since 
they are roots of N y (g,p) = 0, but they are simple poles of ^r(g,p). 

First consider the residues of e pr F(£,p) and £ pT if r (£,p) at the poles ±ia n and ±ib H . The sums of the residues at these 
poles provide the results 


and 


— /> \ "s'' ^ I'fc'ZU-n) ia m r | Nrfj, Jb.r , N rfe H>n) 

m,T) + DA-iaJ e + D Y '(ib„) e + D Y \-ibJ e J 

... (f ja,T . N*(£,— iaj -*lt , N*(gib.) iKr N-i(g—ib n ) 

D7m e + wFi^ e + ww e -'"zvi-i&o e J 


(B29) 

(B30) 


where the subscripts on y and f indicate that they constitute only parts — so-called “dynamic” parts — of the complete solu- 
tions for t/(£,t) and \H£,t). In order to obtain the values of and at p= ±ia„ and ±ib n , it is necessary to use the 
fact that, when p equals one of the roots ira, and ±ib„, one of four equations m 1 =7iir, mi=— nw, mj=imr, or rrh = — wot 

368055 — 50 27 
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is satisfied. The question of which value of p corresponds to each equation need not be answered since the following rela- 
tionships can be derived: 


/iVA =f^\ _ G) nT sm 

. \Pr') mi -nr KDy'J^-nr V^rVmj-lsx KDr'J^-inr 2 ta^+p W+ W 

(N£\ /2VV\ =(W*\ = __©[ 7lV+ (l) pa ] C ° S7l,r? 

wi,. « WA-— 2a^+ v \a^+b/) 


(B31) 


(B32) 


As an example, let p~ia n . Then one of the four possi- 
bilities mi—nTj 77ii =z tltVj rn^—imr^ and m^=—imr must 

N t 

occur. But for any of these possibilities, the ratios y=r-> and 

• u Y 

^7 have the forms given by equations (B31) and (B32). 
Thus, 




0 


mr sin mr£ 


Dy'daJ a* 2 (a n 2 —b/) 


C2\ 2 7l 2 7T 2 


j\cj 0 »‘ 


‘] 


1 cos n-K% 


Dv'iiaJ 




(B33) 


The values of the other necessary ratios can he found in a 
similar manner and equations (B29) and (B30) reduce to 


/c* V 4S mr sin mr£ / 

"COS &*T 

COS b n r\ 

\cj dn— bn y 

K CLn 2 

bn 1 ) 


: ) (B34) 


and 

&,(£> r)=2 g ^ 3 _ b t[(^3- 1) COSa.r-^^j-l) COS & B rJ 

(B35) 

Next, consider the contributions of the pole at p= 0. 
Equations (B24) and (B25) may be expanded in a Laurent 
series about p=0. The resulting equations for F(£,p) and 
^(f,p) may be written 




(B36) 

(B37) 


where 0(p) signifies terms of order p or higher. Thus the 
quantities e^Y^p) and e^fap) have simple poles at p=0 
with the residues 


and 


y*t(z,T) g 2 + 3 




B 


(B38) 


(B39) 


respectively. It can be readily verified that equations (B38) 
and (B39) actually constitute the solution of the problem 
when the beam is considered to be loaded statically. 

Lastly, consider the points p— ± i -/R , which are simple poles 
of ^fap). It can be shown that 


r N * 

hm ^ ; = ■ 


V — ► if 


D 


1 

2 R 


(B40) 


so that the sum of the residues of e^Zrfap) at these poles 
provides the additional dynamic contribution 


&«,(£, f)=— -g cos -JRt 

Summing the contributions of all the poles gives 


(B41) 


(c % sin 7i7rf / cos a*r 

cos £>»r\ 

W« v O, 2 



4 / (£, t ) — cos ^fSr+ 


(B42) 


2 S^ 3 

n-1 A* 


3 7i.7r£ T / 6„ 2 \ 


cos &*r — 


/as! 

\nV m 


1^ cos 6 n rJ 
(B43) 


The responses ?/(£,t) and ^(f,r) are seen to have the same 
form as solutions obtained by Williams' method (ref. 2); 
that is, they are the sum of a static part, y,*(£,T) and &<(£,r), 
and a dynamic correction, y^r) and 4* &,(£,?*), which 

is expanded in a series of the natural mode shapes. The 
y- and ^-components of the natural modes have the shapes 
sin m r£ and cos respectively, and there are seen to bo 
two natural frequencies, a* and b H , associated with each 
integer n so that each integer corresponds to two modes. 
From equations (B2) and (B3), 


M{,r)=H+ 225 n 

77-1 


mV 


(; 


On 


aV 


and 


in nir£ ( b n * 
-bn 1 l\ 

) COS 6,rJ 


0 


cos a*r — 

(B44) 


V (£,t) = COS t/R T— 1+2B Xj C0 » 3 —r- \ (cos OnT—COS b H r ) 

n-1 0 n 

(B46) 
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Six terms (71= 3 ) have been used in the summations of 
equations (B44) and (B45) in computing the quantities 
M(1/2,t ) and V (0,t). These results have been plotted in 
figures 5 (a) and 5 (b) up to t= 8 . 

Simply supported beam with an applied ramp-platform 
end moment. — Let the response to unit step moment at 
£=0, obtained above, be designated -Mi(f,r) and 1 j_(£,t). 
Then the response to an arbitrary applied moment M(0,r) 
may be obtained from Duhamel’s superposition integral as 


r) =M(fi,0)Mi (I, r)+j o Mr (£, T—d)Me(O,ff)d0 ) 
V(t,T)=M( 0, 0)F,& r)+ £v: (^r-&)M e (0,ff)dd j 


(B46) 


where the subscript 6 indicates differentiation with respect 
to that variable. Let the applied moment have the ramp- 
platform time history 


M(0,t ) — t (0 ^ r ^ 1) ) 

M(0,r) = l (0*1) j 

Then equations (B46) become 


and 


^?(£>r) = J^ Mi(t-,T—0)dB 

(O^r^l) 

Mi(t,T—6)de 

(r>l) 

F(f,r)=J o Ti (£,T—6)d3 

(O^r^l) 

y^r)=£v^,r-e)de 

(r>l) 


(B47) 


(B48) 


(B49) 


Substituting Mi _and Vi from equations (B44) and (B45) 
gives for M and V 

(O^T^l) 


( a*? ■ 

A sin 6,r" 

Uv ■ 

l ) K J 


A sin 0 , 7 -— aing,(r— 1) 

) On 

/ a* 8 A sin6,r— sin 6 ,(t— 1) ~1 
Wit*' J K J 


and 


( T >1) 

(B50) 


(Oja 1) 


^ 1 1 or. ^COS TlTrJ/sin o,T sin6„r\ 

j-) 

r(S,T)= T— sin VB (r— 1) J— 1+ 

^ ^ cos tht$ sin g,r— sin a»(T— 1) 

sin 6 »t— sin 6 «(t— 1)1 , ^ , N 

K J (T>1) 

(B51) 

As in the preceding case, computations using terms up to 
71=3 in equations (B50) and (B51) have been made for 
M(1/2 ,t) and F(0 ,t) in the range 0^r^8. The results are 
plotted in figures 6(a) and 6(b). 
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EXACT CLOSED SOLUTIONS FOR UNIFORM BEAMS WITH c t =c 3 


INTRODUCTION 

In this appendix some exact closed solutions of Timo- 
shenko 7 s equations are derived for infinitely long uniform 
beams for which Ci=<h. The results for beams of infinite 
length are then utilized to obtain results for beams of finite 
length. This is done either by simply restricting the infinite- 
beam solution to a time interval in which it coincides with 
the corresponding finite-beam solution or by superposing 
infinite-beam solutions in such a way that the result satisfies 
the boundary conditions for a finite beam. 

EQUATIONS 

If the propagation velocities c x and are taken to. be equal 
and if 2=0, equations (14) reduce to 


COf-b^fr — 0 

— 4X s t? T =0 
Ms+u t -T=0 


Vi—^r2Vr—o)=0 


(Cla) 

(Olb) 


(Old) 


where X=— . This system may be further reduced by 

writing a linear combination of equations (Cla) and '(Clb) 
and a linear combination of (Clc) and (Cld) to obtain 


Zs+2\iK T =0 
Z T -\-2\iK ^ — 2\iZ— 0 


(C2b) 


in terms of the complex variables Z=V+2 and 
K=M+2\iv. 

For a beam initially at rest, equations (C2) may be trans- 
formed to (transform 6, table I) 


f^d-2 XipK — 0 
2Xi/c,+ (p-2Xt)f=0 


in terms of the Laplace transforms 


pa 

k(£,2>) = J o 


K£>2>)=J[ 


(C3a) 

(C3b) 




equation 


e~ VT K{^,r)dr. Eliminating k produces the 


V(J>— 2X£)f=0 


which has the solution 

^,p)=A(p)e~ ( ^^+B(p)e^^ (C4) 

where A and B are governed by the boundary conditions. 
If the beam extends to infinity in the positive f-direetion, 
boundary conditions stipulating that the beam remain 
410 


undisturbed at £= co must be satisfied. Thus, equation (C4) 
reduces to 


&,p)=A(p)e 




Substituting this expression into equation (C3a) produces 

*&p)=^Vp(p- 2 x*>~ w,<p ~™ • (C0) 

The quantity A(jp) may now be determined from the remain- 
ing boundary conditions in conjunction with equation (C6) 
or (C6), or both. 

SPECIFIC PROBLEMS 

Infinite cantilever given a step velocity at the root. — If 
the end £=0 of an infinitely long beam is restrained from 
rotating and is given a unit step velocity at time r=0, the 
boundary conditions are i(0,r) = l and <3(0 ,t) = 0. In terms 
of Z and K these conditions may be written 

I.P.[Z(0,r)]=0 

and 

I.P.[ST(0,r)]=2X 

where I.P. designates the imaginary part of the quantity in 
brackets. The boundary conditions may be transformed to 
(transform 7, table I) 

I-P.lf(0,p)]=0 (C7a) 

and 

I.P.[K(0,p)]=| (07b) 

respectively. 

Because of the nature of these conditions, it will be con- 
venient to proceed in the following manner to determine the 
function A(p). From equation (C5), it may be seen that 
A( 3 >) = S( 0,p); thus condition (07 a) establishes at once that 
A(jp) for this problem is a real quantity. At £=0 equation 
(C6) becomes 

k(°,P)=|^ V?(p-2Xt) 

or, in view of condition (C7b), 

K.P. [ K (0 >P )]+^=4M Vp(p-2Xi) (08) 

where R.P.[x(0,p)] is the real part of *(0,p). Since A(p) and 
p are real, the conjugate equation, which must also hold, is 


R.P.[*(0,p)]- 


2\i A(p) 


Eliminating R.P.[k(0,p)] between equations (C8) and (C9) 
results in 

A(r)=- 2W(V^-V^) (O10) 
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If A($) is substituted into equation (C5), the result may 
be written in the form 

r _ p 5S) 

f (£>P) — — [4A 8 +2Xf (p— Vp 2 +4X g )] ^ === 

which has the inverse transform (transforms 1, 3, 14, and 15, 
table I) 

iJ>Jo(x#=F) ctej 

(r>l) (Oil) 

where Jo and J t are Bessel functions of the first kind. Prom 
equation (Oil) the shear and angular velocity are obtained as 

V (f,r)= — 4X 2 jcos \ t J 0 (xV?=T 2 ) + 

J^sin M J 0 (A#=?) dd^ 

(r>£) (012) 

and 

w(£,r)=-2xjsin XrJoCxVr 1 -? 5 )- 
J £ r cos x^ JoCx# 11 ?) 

(r>?) (013) 

in terms of integrals that apparently defy evaluation by 
analytical methods. 

The time history of the shear 1^(0, r) at the root of the 
beam may be obtained as a special case of equation (C12) ; 
however, it will be noted that L{ V(0 } r) }=f(0,p)=^L(p) 
(where L( F(0 ,r) } denotes the Laplace transform of F(0 ,t)) 
so that the inverse of equation (CIO) gives V(0 ,r) directly. 
Thus the shear at the root is found explicitly in terms of 
tabulated functions as (transform 13, table I) 

T^(0,r) = — 4X 3 [c0S Xr t7o(Xr) +sin Xr Ji(Xr)] (r)>0) (C14) 

Since A(p) may be replaced by i{y(0,r)}, the inverse 
transform of equation (C6) may be written for this problem 
in the form (transforms 1, 2, 3, and 17, table I) 

m,r) i[^^F(Q,r-e- J>(Q,r-e)e^ + 

(t>£) 

or, after integration by parts, 

I^,r)=2XiW 0 (xV^?)- J t T {p(0,r-e)+ 

— a sin X(r-0) JMr-d)}\e M J 0 (x^=¥) dd 
1 — U ) 

(r>£) (Cl 5) 


where V (0,r) is defined by equation (C14). Prom equation 
(C15) the moment and linear velocity may be written as 

M(Z,t)=— 2X sin XtJoCxVt 2 -^)— 2xJ^|^ V(Q,r-ff) cos X0- 

sin X^r-^^Wr-fl)] j Jo(\-y[F=?)dB 

(r>© (C16) 

and 

v(^,t)—COB Xr Jo (xy-r 2 — £ 2 ) — V(0,t — 0) sin X0+ 

sin X(r-fl) Ji[X(r-0)] } J 0 (x#=^)<ft 

(r>9 (017J 

again in terms of integrals requiring evaluation by approxi- 
mate methods. 

As in the case of the shear, it will be more convenient to 
determine the moment M (0,r) at the root by means other 
than a reduction of the general relation, equation (C16). 
Since, for this problem, R.P.[/c(0,p)]=L{5f (0,r)}, equation 
(C8) may be written, after substitution for A(p), in the form 

£{M(0,r)}=i (V£W— p) 

The inverse transform of this equation is (transforms 4 and 
14, table I) 

M(0,r)=2\£ de (t> 0) (0180 

But, from reference 17, 

^^=2XJo(2Xr)-^ [Ji(2Xr)] 

and 

£j 0 (2\6)de= rjj 0 (2Xr) +| [ Jx(2At)H 0 (2At) - 
Jo (2Xr)fli (2 Xt)] j 

where H 0 and Hi are Struve functions. Thus, equation (C18) 
becomes 

M(0,r)=2X{Xr[2-7rff 1 (2Xr)]J o (2Xr)-[l- 1 rXrBo(2XT)]J 1 (2Xr)} 

. (t>0) (019) 

The Struve functions H 0 and Hi are tabulated in reference 
17 in the range 0 g 2 Xt ^ 15.9. . Por larger values of the argu- 
ment it is convenient to use the approximations (ref. 17) 

Ho(2Xt) « Po (2Xt) 4 — r— (2Xr>15.9) 

7TAT 

Hi(2Xt) «Px (2Xt)+| ( 1+4^5) (2Xr>15.9) 

where T 0 and Y x are Bessel functions of the second kind. 
From reference 18, 

Po(2Xr) Jx (2 Xt ) — Yi (2Xr) J 0 (2Xt)=— 

kKt 
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M(0,t)«2x[i-^^ 1 (2Xr>15.9) (C20) 

Computed values of T?(0 ,t) and M(0 ,t), obtained from 
equations (C14), (C19), and (C20) for a beam with X=5, 
have been plotted in figures 4 (a) and 4 (b ), respectively, in 
the range 0 ^r^2. In this range, the root of a finite canti- 
lever beam behaves as if the beam were infinite, since the 
effect of the free end is not felt until the return of the wave 
front at-r=2. 

Infinite simply supported beam with an applied step end 
moment. — If the end £=0 of an infinitely long beam is simply 
supported and is subjected to a unit step bending moment 
at time r=0, the boundary conditions are 5?(0,t) = 1 and 
t>(0,r)=0, or lf(0,r)=l. This condition transforms to 
(transform 7, table I) 

K (0,p)=i (C21) 

Substituting equation (C6) into condition (C21) reveals that 


A(p)=- 


2 Xi 


Vp(p— 2 Xt) 
so that equations (C5) and (C6) become 


r(f,p)=2Xi 


.e 


and 


Vp(P 2Xi) 


( 022 ) 


(C23) 

(024) 


The results of this section are utilized in the next section 
where a finite simply supported beam is considered. 

Finite simply supported beam with an applied step end 
moment. — If the end £=0 of a finite simply supported beam 
is subjected to a unit step bending moment at time r==0, 
the boundary conditions may be written 


M(0,t) = 1 

?(0,r) = M(l,r) = F(l,r) = 0 


(C31) 


or, in terms of K , 


*( 0,0 = 1 ) 
*(1,0 = 0 ) 


(C32) 


Although a direct solution of equations (C2) is possible in 
this case, the response of the finite beam may be obtained 
by a somewhat simpler procedure in which the responses of 
infinite beams to the same disturbance are superposed. 
This procedure is described below. 

A series of semi-infinite simply supported uniform beams 
extending in opposite directions are shown in figure 7. The 
beams have been positioned in the figure so that the origins 
of the space coordinate £ lie on the same vertical line, and 
the segments of greatest interest, 0^£^1, have been out- 
lined with solid lines while the rest of each beam is defined 
by dashed lines. For each beam the origin of the coordinate 
£ is seen to lie at a different position relative to the end whore 
a unit positive moment is suddenly applied at time r=0. 
An infinite number of these beams is assumed to exist. 

The response of the top beam has already been determined 
as (eqs. (C25) and (C26)) 


respectively. The inverse transform of equation (C23) is 
(transforms 3 and 15, table I) 

Zfer)=2Xte^J 0 (xV? =I F) (r>0 (025) 

and the inverse of equation (C24) may be written (transforms 
2, 3, 4, and 16, table I) as 

de (r>0 (026) 


Z($,r)=2X^Jo(xV^ I T 2 ) (r>9 

and 

*(£,T) = ff(£,r) (r>£) 

where 

<?&r)=^— x*J > Jl de 

The response of the second beam may be determined from 
these equations simply by accounting for the changed 


From these equations, the shear, angular velocity, moment, 
and linear velocity are obtained as 


y($ >r )=— 2X sin Xr Jo (X-V? 11 ?), 

(t>9 

(027) 

♦ 

C0 (£>7') = cos Xt«7o (X-^/t 2 — £ 2 ) 

(v>0 

(C28) 

M(£,r)=cos X£ — X£j 

r«. 

(r>9 

(029) 

?($,t)= sin X£ — X^J 

1 { ye 5 -* 2 

(r>© 

(030) 


Thus, V and 5 are determined everywhere in dosed_form in 
terms of tabulated functions, but the relations for M and v 
contain integrals which must be evaluated by approximate 
methods. 


( a — r r t f” — “r 3 

0 12 3 4 



i 

u 


Figure 7. — Semi-infinite uniform beams superposed to obtain the 
response of a finite simply supported beam. 
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direction and shifted origin of the coordinate £. This process gives 

Z(£,t)=2X^Jo [xVt 2 -(2-S) ! J 


and 


m,T)=-G(2-z,T) 


(r>2 -Q 


(r>2-Q 


In order to obtain the response of the third beam from the response of the first, only the shifted origin need be accounted 
for. Thus, for this beam, 

Z(f,r)=2Xw^Jo [xVr s -(2+0 2 ] (t>2+£) 


Similarly, for the fourth beam, 


m,r) = OC. 2+£,r) 

Z(f , r)=2Xte <Xr Jo [xVr 3 -(4-0 s ] 

K(£,t) = — 0(4— £,r) 


( r >2+£) 

(r>4-0 

(r>4-f) 


and so forth, for all the rest. In each case, the region where the response is not specified is a region of zero response. 

Let the responses of all these beams to their respective disturbances be superposed and consider only the response 
of the segment of the resulting composite beam. Since the wave fronts of all the disturbances travel at the 


same velocity ^-=1, the response of this segment may be written 


Z(£,t)=0 

Z(£, T)=2Xie^Jo(xVT s -? 2 ) 

Z(f,r)=2X^{jo(xV^)+Jo[xVr J -(2-O i ]} 

Z(5,r)=2X^{jo(xV : ^)+^o[xVr J -(2-t) 2 ]+Jo[xVr s -(2+f) s ]} 


(O^r^Q' 

(^t^2-|) 

(2-£^t=S2+£) 

(2+^r^4-0 


HC33) 


and 


m,r )= 0 

m,r)=0{k,r) 
m,r)=0^r)-0( ;2-f,r) 
m,r)=0(^r)-0(2-^r) + 0(2-b^r) 


(O^rrgf)' 

(2-^r^2+f) 

(2+f^r^4-f) 


(C34) 


Since 0(O,r)=l, equations (C34) are seen to satisfy the boundary conditions (C32) and the response obtained above 
must be that of a simply supported uniform beam to a unit dimensionless step moment applied at the support {=0. 

The vertical shear V and bending moment M may be obtained from equations (C33) and (C34) as 

VUt) = 0 (Orgr^f) “ 

V(It)=-2\ sin Xr JoCxV^ 1 ! 5 ) (^r^2-Q 

F(£,t)=- 2 X sin Xr{jo(xV^ =I 5 S )+0o[xVT 2 -(2-D a J} (2-^r^2+0 ^° 35 ) 

V(£,t)=-2\ sin Xr{Jo(xV^ 5= y 2 ) + Jo[xVr 1 -(2-?) a ]+J 0 [xVT 2 -(2+Q i j} (2+$^ r^4-£) 


and 


M(£,t) = 0 

M(£,t)=H(£,t)—H(2— £, t ) 

M(f,r) =H(£,t) —H(2—£,t) -\-H(2-\t£,t) 


(Orgr^f)' 

(£^t3S2-£) 

(2-^rg2+{) 

(2+^r^4-$) 


-(C36) 
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.Z7(£,t)=cos X£— X£ f cos 

J C -ye — £ 

The angular and linear velocities S and t? may, of course, 
also be obtained from equations (033) and (C34). _ 

Computations have been made for the quantity T(0,r) 
from equation (C35) with X=5. The results have been 
plotted in figure 5(a) in the range 0^r^8. 

Finite simply supported beam with an applied ramp- 
platform end moment. — Let the applied moment at the end 
£=0 of a simply supported uniform beam have the time 
history 


57(0, r) — T 

57(0, r) = 1 


(O^r^l) ) 

( T >1) ) 


(C37) 


The response of the beam to this disturbance may be ob- 
tained from the response to a unit step by using Duhamel's 
superposition integral. Thus, if the response to a unit step, 
as given by equation (C33), is designated Z i(£,r), the re- 
sponse Z(£,r) to the applied moment (C37) may be written 


ZQzjt) — Zi(£,T—d)d6 
Z(£,r)=J^ Zi(£,r— d)dd 



(C38) 


Substituting from equations (C33) and letting £=0 results in 
the following expressions for the end response Z(0,r): 

Z(0,r)=2XiJ^ T e <x(r -® Jo[X(t— e)]de (0^ 1) 

Z(0,T)=2\i£e ix(T ~ e> J 0 [X(r— d)\dd (l^r^2) 

Z(O ) r)=2XiJ 1 e a(T -»Jo[X(r-e)]d0+ 

4XiJ^ j 0 [xy(r— 0 ) 2 — 4 j dd (2^t^3) 

Z(0,r)=2XiJ o 1 e ^^-» Jo[\(r-d)\ d6+ y (C39) 

■4Xij; C — J 0 [x V(t— 0 ) 2 — 4J de ( 3 ^t^ 4 ) 

Z(0, r)=2\i Jo [X(r- 0)] dB+ 

4Xif>-« J 0 [xV(t — 0) s — 4j dB-\- 
4Xi f r J 0 [xV(r- 0 ) 2 - 16 j de (4^r^5) 


It can be seen that (transforms 3, 4, and 9, table I) 


;{/W°(x*)4 


i d 

~\dp\y v ) 




p-y/p(p—2\i) 

or (transforms 5 and 11, table I) 

e™ Jo (X0) dd— re^ [ Jo (Xt) — i J\ (Xr)] 

Then, Z(0,r) may be reduced in the region 0^r^2 to 
Z(0,T)=2\iJ\™ J 0 (X0) (16 

Z(0,t) =2\ire^[Jo (Xr) - %Jy (Xr)] (O^rrgl) (CftO) 

and 

Z(0,T)=2\i£ e™J o (\0)dB 

=2Xi^ r e* M J°(X5)cte— l e™J 0 QJ)deJ 

Z(0 J T) = 2\ire tXT [Jo(\T) — ft/^Xr)] — 

2Xi(i — l)e**-» { J 0 [X(r- 1 )]-^Ji[X(t- 1)] } 

(l^r^2) (C41) 

From these equations, the shear at £=0 may bo written in 
this range as 

V (0,r) = 2 Xt[co 8 Xr Ji(Xr) — sin Xr Jo(Xr)] (O^r^l) 

V (0,r) = 2Xr[cos Xr Ji(Xr)—sin Xr Jo(Xt)]— 

2X(r— 1) {COS X(r— 1) Ji[X(r- 1)]- 

sin X(r l)t/ 0 [X(r 1)]} (l^r^2) 

In the range r>2, however, the integrals in equations (039) 
apparently cannot be evaluated analytically. If the shear 
F\(0,r) due to a unit step moment has been computed, it 
may be convenient to write, for the shear at the point £=0 
resulting from the disturbance defined by equation (037), 


(042) 


y(0,T)=J Vi(0 ,r — 6)d6 (O^r^l) 

y(°,r)= fy,(0,r-e)d& (r>l) 


(043) 


and then evaluate the integral J V 1 (0,T—d)dd numerically 

to obtain T(0 ,r) in the range r>2. 

The quantity *F(0,r), for a beam with X=5, has been 
computed in the range 0^r^2 from equations (C42) and 
has been obtained in the range 2^r^8 by a numerical 
integration of the exact curve of figure 5 (a) in accordance 
with equation (C43). These results are presented in figure 
6 (a). 
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